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The topological classification of a series of frustration-induced spin-gap phases in the spin-1/2 
ferromagnetic-antiferromagnetic alternating Heisenberg chain with next-nearest-neighbour interaction re¬ 
ported in J. Phys. Soc. Jpn. 82, 064703 (2013) is confirmed using two kinds of entanglement spectra dehned 
by different divisions of the whole chain. For the numerical calculation, the iDMRG method is used. The 
results are consistent with the valence bond solid picture proposed in the previous paper. 


1. Introduction 

Since the discovery of the Haldane phase in spin- 
1 Heisenberg antiferromagnetic chains,^’^ the spin-gap 
phases in the ground states of quantum spin chains 
have been extensively studied in theoretical condensed 
matter physics. Recently, the spin-1 Haldane phase has 
been attracting renewed interest as one of the symmetry- 
protected topological (SPT) phases of matter.It is 
found that this phase is protected by space inversion, 
time reversal, and rr-rotations around two axes.^’'^ From 
this viewpoint, the spin-gap phases of the quantum spin 
chains are classified into SPT and trivial phases in the 
presence of appropriate symmetry. It has also been pro¬ 
posed that these two phases can be distinguished by the 
even-odd parity of the degeneracy of their entanglement 
spectra (ES).^’'* 

In the spin-1/2 frustrated ferromagnetic- 
antiferromagnetic alternating chain with next-nearest- 
neighbour interaction, the present author and coworkers 
found successive phase transitions in which a series of 
trivial and topological spin-gap phases alternate with 
the strength of frustration.® Such successive transitions 
between a series of different spin-gap phases are rather 
unusual and nontrivial. The exact solution is available 
on the ferromagnetic-nonmagnetic phase boundary.®’ 

In the previous paper,® we identified each phase from 
the number of additional S = 1/2 spins necessary to 
compensate the edge spins. This method is, however, 
based on a rather intuitive physical argument. Also, 
there remained ambiguity in the strength of coupling 
between the additional spins and the chain. It is the 
purpose of this work to confirm that these phases are 
actually characterized by the degeneracy of the ES that 
is a well-founded bulk property.®’^ 

This paper is organized as follows. In the next section. 


the model Hamiltonian is presented and the ground-state 
phase diagram obtained in Ref. 8 is reviewed. In sect. 3, 
the entanglement spectrum is introduced and numerical 
results are presented. The last section is devoted to sum¬ 
mary and discussion. 

2. Model and Ground State Phase Diagram 

We consider the S' = 1/2 Heisenberg chains described 
by the Hamiltonian 

L 

'H = ^ {JYS21-1S21 + JAS21S21+1) 

i=i 

2L 

+ ^JlS,S,+2, (1) 

i=l 

where Si is the spin-1/2 operator on the i-th site. In 
this work, we focus on the case Jf, Jl < 0 and Ja > 
0. The lattice structure is shown in Fig. 1. This model 
has the same symmetries as the spin-1 Heisenberg chain, 
if pairs of neighbouring spins are regarded as building 
blocks of the chain. Hence, the same classification of spin- 
gap phases should apply. 

Figure 2 shows the ground-state phase diagram ob¬ 
tained in Ref. 8. Between the Haldane (H) phase and 
the ferromagnetic (F) phase, a series of spin gap phases 
(Is^) are present. In the Is^ phase, the edge spins with 
magnitude S'e appear at the two ends of the open hnite 
chain expressed by the Hamiltonian (1). Each of these 
edge spins is compensated by the addition of 2Se spins 
with magnitude 1/2 on the edge. 

3. Entanglement Spectrum 

To calculate the entanglement spectrum, we divide the 
whole chain into left and right subsystems. The density 
matrices of the right and left subsystems pR and pL are 
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Fig. 1. Lattice structure of the present system. In the calculation 
of the ES, the whole system is divided into left and right subsys¬ 
tems, as indicated by the dotted line (division A) and dash-dotted 
line (division B). 



Fig. 2. Ground state phase diagram. Phase boundaries are taken 
from Ref. 8. The ferromagnetic phase, the Haldane phase, and the 
intermediate spin-gap phases with edge spin Se in the open chain 
are indicated by F, H, and 15 ^, respectively. The dotted line is the 
stability limit of the ferromagnetic phase. The open squares are 
the boundary between different phases determined in Ref. 8. 
The filled circles are the “special points” defined by Dmitriev et 
with exact solutions. The solid curves are guides for the 
eye. The entanglement spectra are calculated along the horizontal 
dotted lines (a)-(d). 


defined by 

PR = Ul |G) (G| (2) 

PL = Tr«|G)(G|, (3) 

where |G) is the ground state of the whole system and 
TrL(R) implies the trace over the left (right) subsystem. 
The eigenvalues Wa and eigenstates |a)R(L) of the density 
matrix pr,(l) satisfy the eigenvalue equations 

PR |a)R = |a)R (4) 

PL |a)L = Wa |q;)l ■ (5) 

The set of eigenvalues {wq.} is common to pr and pl- 
They form an ES of the ground state of the whole system. 
Note that the ES depends on how the whole system is 
divided into subsystems as described below. 


We employ the iDMRG^^’^^ method to calculate the 
ES. In this method, the matrix product wave function 
converges to that of the infinite size fixed point. Hence, 
the effect of the open edge is not reflected in the ES. 
Instead, the degeneracy of the ES results from the virtual 
free spins that appear at the boundary between the left 
and right subsystems. Hence, this degeneracy of the ES 
is a bulk property. 

As we pointed out in Ref. 8, the definitions of the terms 
“topological” and “trivial” are rather arbitrary in the 
present model. In the presence of one of the three symme¬ 
tries that protect the Haldane phase, the ground states 
of the spin chains are classified into those with even and 
odd degeneracies of the ES. The spin-1 Haldane phase, 
which is unambiguously an SPT state, belongs to the for¬ 
mer.^ Our model tends to the spin-1 chain in the limit 
Jf —> —oo, where S 21-1 and S 21 form a spin-1 opera¬ 
tor. In the absence of frustration, the ground state of 
the present model is the Haldane phase if the spin pairs 
connected by Jp bonds are regarded as building blocks 
of the chain. Hence, the division of the spin-1 chain 
at one of the bonds corresponds to the division of the 
chain between S 21 and S 21+1 along line A in Fig. 1 (divi¬ 
sion A) in the present model. Therefore, we identify the 
phase with even ES degeneracy as the SPT phase and 
that with odd ES degeneracy as the trivial phase with 
division A. Thus, to maintain the correspondence with 
the spin-1 chain, we have to adopt division A. 

Nevertheless, one can also define the ES by the division 
of the chain between S' 2 /-i and S 21 along line B in Fig. 
1 (division B). In this case, no virtual free spins appear 
between the left and right subsystems in the Haldane 
phase. In this context, the Haldane phase in our model 
can also be regarded as a trivial spin-gap phase. 

We calculate the ES along the horizontal lines (a)-(d) 
in Fig. 2 for both divisions A and B to obtain deeper in¬ 
sight into the ground state of our model. In most cases, 
the number of the states kept in each subsystem is 240. In 
the neighborhood of the phase boundaries, it is increased 
up to 480 to ensure the convergence. Several large eigen¬ 
values Wa are shown in Fig. 3 for division A and in Fig. 4 
for division B. The degeneracies of the largest eigenvalues 
are also shown in Fig. 3 as gmax for division A and in Fig. 
4 as gmax for division B. In all Ig^ phases, the relations 
gmax = 2S'e +1 and gmax = hold. Physically, this im¬ 
plies that the virtual free spins with magnitude Se appear 
for division A, while those with magnitude S'e — 1/2 ap¬ 
pear for division B at the boundary between the left and 
right subsystems. This confirms that the /g, phases with 
integer Se and half-odd integer Se are topologically dis¬ 
tinct. As for the valence bond structure, this also means 
that the number of valence bonds is 2Se across division 
A, while it is 2Se — 1 across division B. This is consistent 
with the valence bond solid picture proposed in Ref. 8. 

Unfortunately, the convergence of the iDMRG is not 
good near the phase boundary. Hence, it is difficult to 
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Fig. 3. Jp-dependences of several large eigenvalues Wa. of the 
density matrices Pr(l) ^-^d the degeneracy of the largest eigenvalue 
pmax with division A for (a) Jp = —1.75, (b) —2, (c) —2.25, and 
(d) —2.5 with Ja = 1. The vertical dotted lines are the phase 
boundaries determined in Ref. 8. 


determine the phase boundary using only the iDMRG 
results. Although some of the data for the phase bound¬ 
aries taken from Ref. 8 are slightly inconsistent with the 
present iDMRG data, these would be attributed to the 
ambiguity in the size-extrapolation procedures in Ref 8. 
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Fig. 4. Jp-dependences of several large eigenvalues Wa. of the 
density matrices Pr(l) ^rid the degeneracy of the largest eigenvalue 
pmax with division B for (a) Jp = —1.75, (b) —2, (c) —2.25, and 
(d) —2.5 with Ja = 1. The vertical dotted lines are the phase 
boundaries determined in Ref. 8. 


4. Summary and Discussion 

We employed the iDMRG method to calculate the ES 
of the spin-1/2 ferromagnetic-antiferromagnetic alternat¬ 
ing Heisenberg chain with next-nearest-neighbour inter- 
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action for divisions A and B. We have confirmed the 
presence of successive frustration-induced phase transi¬ 
tions between a series of topological and trivial spin-gap 
phases predicted in Ref. 8 using the characterization of 
these phases in terms of the ES. The valence bond solid 
picture of each phase is confirmed. 

In the frustrated quantum spin chains, various kinds of 
exotic phases induced by the interplay of quantum fluc¬ 
tuation, frustration, and ferromagnetic correlation are 
predicted near the frustration-induced transition point 
between the ferromagnetic and nonmagnetic phases. For 
example, the spin-nematic phase^^“^® and partial fer- 
rimagnetic phases with incommensurate spin correla¬ 
tion^^’ been extensively discussed. The pres¬ 

ence of successive transitions between SPT and trivial 
spin-gap phases is also one of the exotic phenomena in 
frustrated quantum spin chains. Recently, these succes¬ 
sive transitions have also been confirmed by the mapping 
onto the nonlinear a model.The possibility of sim¬ 
ilar successive transitions has been suggested in other 
frustrated spin chains,although the results were not 
conclusive owing to the limited system size. Further theo¬ 
retical and experimental searches for these exotic phases 
would be promising subjects in the field of frustrated 
quantum magnetism. 

The author thanks S. C. Furuya for enlightening com¬ 
ments and discussion on the earlier version of this paper. 
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ter, Institute for Solid State Physics, University of Tokyo, 
and Yukawa Institute Computer Facility in Kyoto Uni¬ 
versity. This work is supported by a Grant-in-Aid for 
Scientific Research (C) (25400389) from the Japan Soci¬ 
ety for the Promotion of Science. 
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